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ABSTRACT 

We present a theoretical investigation of the effect of multiple ionisation sources in H n 
regions on the total elemental abundances derived from the analysis of collisionally excited 
emission lines. We focus on empirical methods based on direct temperature measurements that 
are commonly employed in cases when the temperature of the nebular gas can be determined 
from the ratio of nebular to auroral lines of (e.g.) doubly ionised oxygen. We find that direct 
temperature methods that employ a two-temperature zone approach (DT2T methods) are very 
robust against the spatial distribution of sources. Errors smaller than 0. 15 dex are estimated for 
regions where the metallicity is twice solar and errors below 0.05 dex for solar metallicities 
and below. The biases introduced by the spatial distribution of the ionisation sources are thus 
much smaller for DT2T methods than for strong line methods, previously investigated by 
Ercolano, Bastian & Stasihska. Our findings are in agreement with the recent study of H n 
regions in NGC 300 by Bresolin et al. 

Key words: galaxies:abundances;galaxies:ISM;ISM:HII regions;ISM:abundances 



> 

m 
in 
m 

ON 

o 

&\ 
O 



1 INTRODUCTION 

The analysis of emission lines from H n regions powered by OB 
stars is often the only means available for the determination of gas 
abundances both in our Galaxy and others. Accurate abundance de- 
terminations are crucial to derive metallicity gradients, which pro- 
vide a key observational constraint to chemical evolution models 
of galaxies. The spectra of H n regions in the optical and infrared 
are dominated by collisionally excited lines (CELs) of singly and 
doubly ionised ions of some of the more abundant heavy elements 
(e.g. oxygen, carbon etc.). While the metallicity diagnostic power 
of CELs from H n regions is widely recognised and used for study- 
ing the chemical composition of galaxies, the uncertainties inherent 
to the empirical methods employed are often overlooked, poten- 
tially leading to somewhat over-optimistic error estimates. 

Ionic abundances can be obtained from emission lines via 
the solution of the statistical equilibrium equation under the as- 
sumption of a gas temperature and density. Both gas tempera- 
ture and density can be empirically derived from the observa- 
tions using the ratio of diagnostic emission lines (e.g. Osterbrock 
and Ferland, 2006). A widely used temperature diagnostic re- 
lies on the ratio of nebular to auroral [OIII] lines, in particular 
([Om]^5007+[Om]i4959)/[Om]i4363. This method is often re- 
ferred to as the 'direct temperature' method. For faint (distant) or 
metal-rich regions, however, the [OIII]/14363 line is seldom de- 
tected, rendering impossible the direct determination of tempera- 
ture. In these cases one is often forced to use 'strong-line methods', 
which rely on ratios of some of the strongest CELs calibrated with 



one-dimensional photoionisation models. In a previous paper (Er- 
colano, Bastian & Stasihska, 2007, EBS07) we showed that com- 
monly employed strong-line methods, may introduce a systematic- 
bias of typically 0.1-0.3 dex, but up to 1 dex on the logarithmic 
oxygen abundance for regions where the gas and the ionising stars 
are intermixed. 

Direct-temperature methods should be more reliable than 
strong line methods; however even these are not entirely immune 
by errors introduced by calibrations from spherically symmetric 
models. While the error on the determination of a single ionic abun- 
dance from emission lines via the solution of the statistical equilib- 
rium matrix is simply determined by the accuracy of the gas tem- 
perature and (to a lesser extent) density estimates, not all the ionic 
stages of a given element produce emission lines in the observed 
wavelength range. For this reason the determination of total ele- 
mental abundances from ionic abundances relies on a correction for 
the unseen stages of ionisation, known as the ionisation correction 
factor (ICF). A potential bias is therefore introduced by the ICF 
schemes that are themselves calibrated via one-dimensional pho- 
toionisation modeling (e.g. Kingsburgh & Barlow 1994, Peimbert 
et al., 1992), with the implied assumption of a geometry comprising 
of a spherical region with a single central location for all ionisation 
sources. 

The aim of this paper is to estimate whether the use of ICF 
schemes to derive total elemental abundances of H n regions from 
direct-temperature methods introduces a bias in the case of regions 
ionised by multiple stars intermixed with the nebular gas. Via the 
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analysis of theoretical H n region spectra obtained with the mo- 
cassin code (Ercolano et al. 2003, 2005, 2008) in the set-up used 
by EBS07, we show that, unlike strong line methods, direct tem- 
perature methods using a two-temperature approach are very reli- 
able also in cases when the stars are fully distributed within the 
nebular gas. These results are in line with the recent findings of 
Bresolin et al (2009, B09) who found good agreement between 
CEL abundances obtained by direct temperature methods with the 
abundances derived from absorption line analysis of stellar photo- 
spheres in a sample of H n regions in the spiral galaxy NGC 300. 

The paper is organised as follows. Section 2 summarises the 
model setup and input parameters. Section 3 contains a brief de- 
scription of the methods employed. Our results are given in Section 
4, while Section 5 is dedicated to a discussion and conclusions. 



2 THEORETICAL EMISSION LINE SPECTRA: MODEL 
SETUP AND INPUT PARAMETERS 

We have used the theoretical nebular spectra of EBS07, which 
were obtained using the three-dimensional photoionisation code 
mocassin (Ercolano et al. 2003, 2005, 2008). This code uses a 
Monte Carlo approach to the transfer of radiation and can easily 
deal with multiple ionisation sources arbitrarily distributed within 
the simulation region. The atomic database included opacity data 
from Verner et al. (1993) and Verner & Yakovlev (1995), energy 
levels, collision strengths and transition probabilities from Version 
5.2 of the CHIANTI database (Landi et al. 2006, and references 
therein) and the improved hydrogen and helium free-bound contin- 
uous emission data of Ercolano & Storey (2006). The model setup 
and input parameters are described in EBS07, here we summarise 
briefly the main points, but refer the interested reader to EBS07 for 
further details. 

The regions are assumed to be spherical and consisting of 
homogeneous gas with number density N# = 100 cirr 3 . The to- 
tal number of ionising photons is constant for all models and is 
Q H o = 2.8 x 10 50 s . We consider models of five different metal- 
licities ( Z/Z Q = 0.05, 0.2, 0.4, 1.0 and 2.0). The solar abundance 
model assumes the values of Grevesse and Sauval (1998) with the 
exception of C, N and O abundances which are taken from Allende 
Prieto et al. (2002), Holweger (2001) and Allende Prieto, Lambert 
and Asplund (2001), respectively. The higher and lower metallicity 
cases were obtained from the solar abundances by scaling using the 
empirical abundance trends observed in H n regions by Izotov et al 
(2006). The gas is ionised by 240 sources belonging to two pop- 
ulations, a hot (M, = 56 M Q ) and cool (M« = 37 M Q ) population, 
each population, as a whole, emits equal quantities of H-ionising 
photons. The stellar spectra were computed with the starburst99 
spectral synthesis code (Leitherer et al. 1999) with the up-to-date 
non-LTE stellar atmospheres implemented by Smith, Norris and 
Crowther (2002) using single isochrones for the appropriate stellar 
masses. The stars were distrubuted as follows: (i) centrally con- 
centrated at the centre of the spherical region - C-models; (ii) dis- 
tributed in the half- volume of the spherical region - H-models; (iii) 
distributed in the full spherical volume - F-models. In the C-models 
all stars share the same location at the origin of the Cartesian axes. 
In the H- and F-models the stars are distributed stochastically such 
as to obtain a statistically homogeneous 3D distribution of sources 
in the half or full spherical region, respectively. The Stroemgren 
sphere of stars in the F-models seldom overlap, while those of the 
C-model completely overlap, with the H-model representing the in- 
termediate case. 



In Table 1 we list the subset of the emission lines from our 
theoretical spectra that were used for the analysis described in the 
following section. 



3 THE DIRECT TEMPERATURE METHOD 

The aim of this study is to test whether the abundances determined 
from the emission line spectrum of an H n region ionised by multi- 
ple stars intermixed with the gas are different from the abundances 
determined by a region with exactly the same physical character- 
istics but with all ionisation sources concentrated in the centre. 
EBS07 indeed demonstrated that large biases are introduced by this 
effect when strong-line methods are used. The main reason for such 
differences was due to a decrease of the 'effective ionisation pa- 
rameter' of the gas when the stars were fully distributed within the 
medium compared to when the same stars were all concentrated at 
the centre of the nebula. This affected the temperature structure of 
the nebula significantly enough to produce large errors in the de- 
rived abundances. 

Ionic abundances from direct-temperature methods should be 
immune from this error as long as the temperature gradients within 
a given ionic phase are not too large (see e.g. Stasinska, 1980 and 
Kingdon & Ferland, 1995). However total elemental abundances 
can only be obtained by applying an ICF scheme to correct for 
the unseen ionisation stages. The question therefore remains as to 
what is the effect of the geometrical distribution of the stars on the 
ICFs, which rely on theoretical calibrations via one-dimensional 
photoionisation models (e.g. Kingsburgh & Barlow, 1994, hereafter 
KB94; Peimbert et al., 1992, hereafter PG92). 

To answer this question we took a subset of lines typically 
observed in extragalactic H n regions, including the important tem- 
perature diagnostic lines of [O in] at 4363A and [N n] at 5754A, 
from the line spectra produced by the photoionisation models de- 
scribed in Section 2. We then used these model spectra to derive 
chemical abundances via the direct temperature method. The lines 
used for the analysis are listed in Table 1 . 

We considered the nebula as being composed of two separate 
zones of low and high ionisation: first, we assumed a temperature of 
10 000 K to obtain electron densities from the [O n] ^3726//13729 
and [S n] /16717//16731 line ratios. Then, the average of these two 
electron densities was used to derive a temperature from the [N n] 
(/16548+6584)//i5754 line ratio. This temperature was then used 
to recalculate the densities, and the temperature recalculated once 
more using the resulting density. The abundances of singly-ionised 
species were derived using this temperature and density. Then, we 
used the same iterative approach, but using the [CI in] 5517/5537 
and [Ar iv] 47 1 1/4740 line ratios as density diagnostics, and the 
[O m] 4959+5007/4363 ratio as a temperature diagnostic. Abun- 
dances of doubly and more highly ionised species were derived us- 
ing this temperature and density. 

Total abundances relative to hydrogen were calculated for He, 
C, N, O, Ne, S and Ar, using the two commonly used ICF schemes 
of KB94 and PTPR92. For several atoms these schemes use the 
same correction; they differ for helium (KB94 does not correct for 
neutral helium while PTPR92 does), argon and sulfur. 
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Table 1. Model emission line fluxes, relative to H/J=1()0. 



A (A) Ion Model 









0.05c 


0.05f 


0.05h 


0.20c 


0.20f 


0.20h 


0.40c 


0.40f 


0.40h 


1.00c 


l.OOf 


l.OOh 


2.00c 


2.00f 


2.00h 


3726. 


03 


[On] 


12.38 


27.70 


14.07 


40.05 


87.69 


35.55 


35.20 


95.33 


32.05 


18.55 


48.28 


17.57 


63.89 


62.60 


55.92 


3728. 


.82 


[On] 


16.51 


37.10 


18.78 


53.18 


116.7 


47.19 


46.62 


126.5 


42.47 


24.34 


63.43 


23.04 


82.17 


80.68 


71.92 


3868. 


.75 


[Neii] 


22.93 


11.58 


21.12 


43.07 


26.43 


44.01 


31.03 


21.10 


31.42 


2.682 


2.433 


2.410 


0.034 


0.149 


0.040 


3967. 


.46 


[Nen] 


6.910 


3.491 


6.363 


12.97 


7.965 


13.25 


9.349 


6.359 


9.468 


0.808 


0.733 


0.726 


0.010 


0.045 


0.012 


4068. 


.60 


[Sn] 


0.471 


1.828 


0.634 


1.426 


5.011 


1.341 


1.590 


6.435 


1.547 


1.020 


4.000 


0.950 


0.351 


1.072 


0.331 


4076. 


.35 


[Sn] 


0.163 


0.632 


0.219 


0.493 


1.734 


0.464 


0.550 


2.227 


0.535 


0.354 


1.390 


0.329 


0.122 


0.374 


0.115 


4363. 


.21 


[O m] 


9.693 


3.995 


8.741 


8.649 


4.593 


9.138 


2.993 


1.747 


3.037 


0.070 


0.045 


0.053 


0.002 


0.005 


0.002 


4471. 


.50 


He i 


4.150 


4.186 


4.143 


4.328 


4.241 


4.305 


4.469 


4.295 


4.460 


4.723 


4.476 


4.723 


4.293 


4.356 


4.345 


4685. 


.68 


Hen 


0.008 


0.011 


0.010 


0.013 


0.013 


0.013 


0.007 


0.010 


0.009 


0.006 


0.009 


0.006 


0.003 


0.001 


0.002 


4711. 


.37 


[Ar iv] 


1.385 


0.491 


1.235 


2.319 


0.898 


2.451 


1.765 


0.710 


1.890 


0.081 


0.047 


0.082 


0.000 


0.002 


0.001 


4740. 


.17 


[Ar iv] 


1.069 


0.377 


0.952 


1.737 


0.672 


1.838 


1.295 


0.520 


1.386 


0.059 


0.034 


0.060 


0.000 


0.001 


0.000 


4861. 


.33 


Hi 


100.0 


100.0 


100.0 


100.0 


100.0 


100.0 


100.0 


100.0 


100.0 


100.0 


100.0 


100.0 


100.0 


100.0 


100.0 


4958. 


91 


[O m] 


110.1 


52.62 


101.5 


200.4 


116.4 


207.7 


161.3 


102.3 


165.2 


19.47 


16.55 


17.97 


1.854 


4.190 


1.844 


5006. 


.84 


[O in] 


328.6 


157.0 


303.0 


598.0 


347.4 


619.8 


481.5 


305.2 


493.0 


58.11 


49.40 


53.64 


5.531 


12.50 


5.504 


5754. 


.60 


[Nii] 


0.026 


0.059 


0.030 


0.076 


0.165 


0.068 


0.097 


0.258 


0.089 


0.075 


0.193 


0.071 


0.237 


0.266 


0.215 


5875. 


.66 


He i 


10.87 


11.16 


10.87 


11.57 


11.42 


11.50 


12.27 


11.80 


12.23 


14.32 


13.22 


14.26 


14.12 


13.42 


14.15 


6312. 


.10 


[S m] 


0.923 


1.045 


0.982 


2.261 


2.867 


2.142 


1.868 


2.573 


1.686 


0.458 


0.743 


0.410 


0.368 


0.309 


0.320 


6548. 


.10 


[Nii] 


0.321 


0.943 


0.396 


1.114 


2.919 


1.017 


1.895 


5.979 


1.784 


2.831 


9.155 


2.741 


16.14 


21.44 


15.15 


6562. 


.77 


Hi 


276.8 


280.8 


277.3 


281.2 


283.0 


281.1 


286.8 


287.1 


286.7 


306.5 


301.3 


305.7 


314.4 


306.6 


313.7 


6583. 


.50 


[Nii] 


0.981 


2.880 


1.209 


3.402 


8.916 


3.107 


5.787 


18.26 


5.449 


8.647 


27.96 


8.371 


49.31 


65.50 


46.28 


6678. 


16 


He i 


2.983 


3.115 


2.993 


3.238 


3.214 


3.218 


3.486 


3.354 


3.476 


4.080 


3.784 


4.069 


3.864 


3.792 


3.905 


6716. 


.44 


[Sn] 


4.194 


17.80 


5.839 


12.84 


50.82 


12.43 


16.03 


70.59 


15.77 


13.80 


61.42 


13.08 


6.061 


21.11 


5.783 


6730. 


.82 


[Sn] 


3.100 


13.11 


4.314 


9.524 


37.39 


9.199 


11.83 


51.96 


11.64 


10.15 


45.06 


9.619 


4.447 


15.50 


4.242 


7136. 


.80 


[Ar m] 


0.808 


1.193 


0.875 


2.969 


4.414 


2.833 


3.662 


5.865 


3.412 


2.342 


4.033 


2.205 


3.014 


4.011 


2.592 


7319 




[On] 


0.544 


1.025 


0.597 


1.479 


2.856 


1.290 


1.015 


2.393 


0.902 


0.311 


0.676 


0.291 


0.644 


0.568 


0.551 


7330 




[On] 


0.443 


0.836 


0.487 


1.205 


2.329 


1.052 


0.828 


1.953 


0.736 


0.254 


0.552 


0.238 


0.526 


0.464 


0.450 


7751. 


.43 


[Ar hi] 


0.194 


0.286 


0.210 


0.711 


1.058 


0.679 


0.877 


1.405 


0.818 


0.561 


0.966 


0.528 


0.722 


0.961 


0.621 
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Figure 1. Error on the logarithmic abundances calculated via the DT2T 
method caused by the spatial distribution of ionisation sources. The black 
solid line shows results obtained via the KB94 ICFs, while the red dashed 
line shows those obtained via the PTR92 ICFs. 



We also considered the case where a low-ionisation tempera- 
ture diagnostic is not available. In this situation, one can either use 
the [O in] temperature for all ions, or estimate the low ionisation 
temperature using relations such as that found by Pilyugin et al. 
(2006). We confirm previous findings that applying a single temper- 
ature diagnostics introduces considerable errors into the abundance 
determinations. This case will not be further discussed. 



4 RESULTS 

Our analysis indicates that the errors on the abundances intro- 
duced by the geometrical distribution of the ionising stars are much 
smaller when direct temperature methods are used instead of strong 
line methods. We stress here that our aim is not to assess the valid- 
ity of the ICFs but to assess the effect of the distribution of the 
ionising sources on the derived abundances. Therefore, we are not 
comparing the derived abundances to the 'right answer' (i.e. the in- 
put abundances for our models), but rather comparing the results 
from model nebulae with distributed ionising sources to those with 
a centrally concentrated source. 

Figure 1 shows the logarithmic error, E F _ C , on the abundances 
of various elements due to the spatial distribution of the ionisation 
sources. The black solid line shows results using KB94 ICFs and 
the red dashed lines shows the results for the ICF scheme of PTR92. 
The errors are due partly to temperature effects (e.g. a steep tem- 
perature profile or differences in the mean temperatures of various 
ions) and partly to the ICFs which may be more or less sensitive 
to changes in the effective ionisation parameter brought about by a 
different spatial distribution of sources (see discussion in EBS07). 
Here we are mostly interested in the latter, and in order to isolate 
this effect we calculate a correction to the error due to the ICFs, 
A(E F _ C ), by comparing the C and F-model theoretical and empiri- 
cal ICFs according to the following: 

I ICF" ICFl\ 

AC*«) = toKo^/j^) (1) 
where ICF", ICF" are the observational ICFs for the F- and C- 



models, respectively, and ICF T F , ICF^, are their theoretical coun- 
terparts calculated by the photoionisation model. 

The logarithmic errors and the ICF corrections are sum- 
marised in Table 2. In the following we will discuss the sources 
of error in more detail for each element. 



4.1 Helium 

KB94 do not include an ICF correction for the unseen neutral he- 
lium; the very small Ejr-c for KB94 shown in Figure 1 is therefore 
only due to the lack of correction and the C- and F-models hav- 
ing different amounts of neutral helium. PTR92 on the other hand 
correct for neutral helium and it is indeed the ICF correction em- 
ployed that at low metallicities is sensitive to changes in the effec- 
tive ionisation parameter. The A{E F _ c ) corrections given in Table 2 
drastically reduce the E F -c values for PTR92. 



4.2 Oxygen 

The abundance of oxygen derived by the direct temperature method 
with a two temperature description of the medium (DT2T) is 
not very sensitive to the geometrical distribution of the ionisation 
sources (i.e. to the effective ionisation parameter) for the range of 
metallicities discussed here. The largest errors occur for metal-rich 
regions (twice solar) and are always below 0.15 dex. We note that 
both the empirical and theoretical ICFs for oxygen are roughly 
unity for all models considered here. The small error at higher 
metallicities is due to the steepening of the temperature profile 
which is more accentuated for the C-models (see discussion in 
EBS07 and Stasihska 1980). This causes the oxygen abundances 
to be underestimated in the C-models more than in the F-models 
producing the error observed. 



4.3 Nitrogen 

Nitrogen also shows a similar behaviour with a slightly larger error 
(0.21 dex) for metal-rich regions. The nitrogen ICFs are not unity 
and some of the error shown in Figure 1 are indeed due to a dif- 
ferent response of the ICFs to the change in the effective ionisation 
potential. The ICF correction is A(E F _ C ) = - 0.14 for nitrogen 
in the Z/Z Q = 2 case, which brings the nitrogen E F _ c to roughly 
0.1 dex. The remainder of the error can again be ascribed to the 
steepening of the temperature profile as discussed above. 



4.4 Neon 

The situation for neon is more complicated. A(E F _ C ) at Z/Z Q = 2 
is — 1-0.3 dex, which actually increases the magnitude of the error. 
The large discrepancy (E F -c) ~ 0.7!) is due to the displacement of 
the Ne 2+ and 2+ regions which results in a significant difference 
between T e (Ne 2+ ) and T c (0 2+ ) at high Z's. Using T e (0 2+ ) in the 
determination of Ne 2+ abundance then results in the large error. 
Indeed we find that at Z/Z Q = 2, T e (Ne 2+ )/T f (0 2+ ) = 0.85 for the 
C-model and 0.93 for the F-model, the difference has a significant 
impact on abundances derived via CELs due to their exponential 
dependance on temperature. We calculate the correction due to the 
difference in the two temperature regions to be -0.61 dex, which 
brings the total E F -c to roughly 0.1 dex, which is comparable to 
what we found for oxygen and nitrogen. 
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Table 2. Logarithmic errors and ICF corrections 



Helium Oxygen Nitrogen 

Z/Z G £>_ c A(£f_ c ) Z/Z Q £>_ c A(£ F _ C ) Z/Z Q £ f _ c A(£ f _ c ) 

KB94 PTR92 KB94 PTR92 KB94 PTR92 KB94 PTR92 KB94 PTR92 KB94 PTR92 



0.05 


0.01 


0.22 


1.0 


-0.14 


0.05 


0.03 


0.03 


1.0 


1.0 


0.05 


0.07 


0.07 


-0.03 


-0.03 


0.2 


-0.01 


0.11 


1.0 


-0.06 


0.2 


0.03 


0.03 


1.0 


1.0 


0.2 


0.03 


0.03 


+0.08 


+0.08 


0.4 


-0.01 


0.08 


1.0 


-0.03 


0.4 


0.02 


0.02 


1.0 


1.0 


0.4 


0.04 


0.04 


+0.03 


+0.03 


1.0 


-0.03 


0.03 


1.0 


-0.007 


1.0 


0.06 


0.06 


1.0 


1.0 


1.0 


0.08 


0.08 


-0.03 


-0.03 


2.0 


-0.02 


0.02 


1.0 


-0.047 


2.0 


0.13 


0.13 


1.0 


1.0 


2.0 


0.21 


0.21 


-0.14 


-0.14 



Neon Argon Sulfur 

Z/Z £>_ c MEf-c) Z/Z q £>_ c A(£ f _ c ) Z/Z q £ f _ c A(E f _ c ) 

KB94 PTR92 KB94 PTR92 KB94 PTR92 KB94 PTR92 KB94 PTR92 KB94 PTR92 



0.05 


0.06 


0.06 


-0.07 


-0.07 


0.2 


0.06 


0.06 


-0.04 


-0.04 


0.4 


0.06 


0.06 


-0.006 


-0.006 


1.0 


0.05 


0.05 


+0.10 


+0.10 


2.0 


0.39 


0.39 


+0.41 


+0.41 



0.05 


0.15 


-0.11 


-0.07 


1.0 


0.2 


0.20 


-0.03 


-0.04 


1.0 


0.4 


0.23 


0.04 


-0.004 


1.0 


1.0 


0.29 
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4.5 Argon 

PTR92 do not employ a correction for Ar + , and the small errors 
shown are due to the change in the ionisation structure of F- and 
C-models. KB94 do include an ICF for Ar, however this is quite 
sensitive to changes to the effective ionisation parameter. The max- 
imum errors for KB 94 are of 0.3 dex at solar metallicities against a 
maximum error of 0.1 dex obtained by PTR92 at Z/Z Q = 0.05. 



4.6 Sulfur 

Sulfur presents larger problems at low metallicities for both meth- 
ods with maximum errors of 0.4 dex at Z/Z G = 0.05 for PTR92 
and 0.27 dex for KB94 at the same metallicity. Most of the error 
here can be ascribed to the different responses of the ICFs to the 
change of the effective ionisation parameter. Indeed applying the 
A(E F -c) corrections given in Table 2 brings E F -c to values smaller 
than 0.1 dex both for KB94 and PTR92. 

The results obtained for the H-models lie in between those of 
the F and C models and are therefore not included in this discus- 
sion. 



5 DISCUSSION & CONCLUSIONS 

The main conclusion of this short paper is that abundance deter- 
minations from collisionally excited emission lines (CELs) of H n 
regions via direct temperature methods that use a two-temperature 
description of the ionised region are very robust and not affected 
significantly by the spatial distribution of ionisation sources. In- 
deed the maximum errors on the oxygen abundance derived with 
the DT2T method are still lower than 0.15 dex at Z = 2 and below 
0.05 dex at solar metallicity and below. 

As a comparison we found that the strong-line-methods anal- 
ysed in EBS07 gave much larger errors as shown in Figure 2. The 
figure shows the discrepancy in the empirically determined oxygen 
abundance for the F- (fully distributed) and C- (centrally concen- 
trated) models. The thick black line indicates the results for the 
DT2T, while the thinner lines are for the results of the various 
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Figure 2. Errors on the logarithmic abundances calculated via the DT2T 
method caused by the spatial distribution of ionisation sources. The black 
solid line shows results obtained via the DT2T method. The thinner lines 
show results obtained from strong line methods, namely as follows. Red 
dashed: 023 (Pilyugin, 2000, 2001b); green dotted: 03N2 (Stasihska 2006); 
blue dash dotted: N2 (Pettini & Pagel 2004); magenta dash double dot: S23 
(Perez-Montero & Diaz 2005); cyan long dash: S303 (Stasiriska 2006); 
black thin solid: Ar303 (Stasihska 2006). The yellow shaded region in- 
dicates the metallicity range of the H n regions in NGC 300 analysed by 
B09. 



strong line methods considered by EBS07. A detailed legend is pro- 
vided in the caption to Figure 2. The small errors associated with 
the DT2T method compared to the large errors of the strong line 
methods are in agreement with the recent results of Bresolin et al. 
(2009, B09) who found a systematic difference between the oxy- 
gen abundance calculated by strong line methods and direct tem- 
perature methods for a sample of H n regions in NGC 300. The 
yellow section in Figure 2 shows the metallicity range of the H n re- 
gions analysed by B09. B09 compared the metallicities obtained by 
a direct temperature and a strong-line analysis of the emission line 
spectra of H n regions in NGC 300 to those obtained from B and A 
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Figure 3. Errors on the logarithmic abundances calculated via the DT2T 
method caused by the spatial distribution of ionisation sources. The black 
solid line shows results obtained via the DT2T method. The thinner lines 
show results obtained from the strong line methods analysed by B09. See 
text for details. 



only one model data point actually falls in that range. In view of 
these shortcomings of the models we can at present only suggest 
that the spatial distribution of ionising sources is the cause of the 
metallicity bias that afflicts strong line measurements and postpone 
firmer statements to a future work where the parameter range is 
better suited to match those particular observations. 

We finally note that we have not included a discussion of the 
well known abundance discrepancy between CELs and recombi- 
nation lines (RLs). A number of possible causes has been identi- 
fied in the literature including temperature fluctuations (Peimbert, 
1967), hydrogen-deficient, metal-rich inclusions (Liu et al., 2000; 
Stasinska et al. 2007) and X-ray irradiated quasi-neutral clumps 
(Ercolano, 2009). The jury is still out however as to which of the 
above effects or a combination thereof is to blame for the discrep- 
ancy. Until the latter problem is resolved all abundances determined 
via nebular emission lines carry a potential error. The excellent 
agreement between the results obtained by B09 from direct temper- 
ature analysis of CELs and those from the supergiants in NGC 300, 
however, indicates that in this Galaxy temperature fluctuations and 
X-ray irradiated quasineutral clumps, if at all present, must be play- 
ing a minor role. 



supergiants in the same galaxy. They found excellent agreement be- 
tween the results from supergiants and direct temperature analysis 
of the H ii regions, while noticing that a systematic bias affected the 
results from some popular calibration of strong lines. The calibra- 
tions used in B09 were not considered by EBS07 and included (i) 
the R233 ratio from the theoretical calibration of Mc Gaugh (1991, 
M91) using the analytical prescriptions of Kuzio de Naray et al. 
(2004) and Tremonti et al. (2004, their Eq. 1, T04), (ii) the theo- 
retical prediction for the [N H]A6583/[0 II]^3727 ratio by Kewley 
& Dopita (2002, KD02), and (iii) N2 = log([N II]i6583/Ha), cali- 
brated empirically by Pettini & Pagel (2004, PP04). 

In order to estimate whether the spatial distribution of stars 
may be playing a role in producing the bias observed by B09 and 
predicted by ECD07, we have used the emission line spectra in Ta- 
ble 1 to compute the oxygen abundances given by the M91, T04, 
KD02 and PP04 calibrations listed above and compared it to the 
DT2T results in Figure 3. The metallicity range of the B09 sample 
is again highlighted by the yellow section. The errors of the strong 
line methods are comparable to those lamented in the B09 paper, 
however a detailed comparison with the observation is premature at 
this point. The main problem is that the set of models run by EBS07 
comprised a very idealised ionising source population which was 
designed to highlight eventual temperature fluctuations that may be 
introduced by the distribution of stars with spectra of very different 
hardness (see also Section|2j, which turns out to be equivalent to a 
much harder 'effective' spectrum than that inferred by B09 for the 
H ii regions in NGC 300. The parameter r\ = (0 + /0 ++ )/(S + /S ++ ), 
(Vilchez & Pagel, 1988) was introduced as a measure of the hard- 
ness of the ionising field, with larger numbers corresponding to a 
softer spectrum. B09 find an average log(^) parameter of roughly 
0.7, while we find values ranging between -0.2 and 0.2, indicating a 
significantly harder spectrum than that of the B09 H n regions. An- 
other problem is the fact that EBS07 explored a wide metallicity 
range and as a consequence the narrow metallicity range of the H n 
regions in NGC 300 is very sparsely sampled, as shown in Figure 3 

1 R23 = ([O H]-13727 + [O III]i37274959,5007)/HyS, 
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